Non-monotonic dependence of the rupture force in polymer chains on their lengths 
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We consider the rupture dynamics of a homopolymer chain pulled at one end at a constant loading 
rate. Our model of the breakable polymer is related to the Rouse chain, with the only difference that 
the interaction between the monomers is described by the Morse potential instead of the harmonic 
one, and thus allows for mechanical failure. We show that in the experimentally relevant domain 
of parameters the dependence of the most probable rupture force on the chain length may be non- 
monotonic, so that the medium-length chains break easier than the short and the long ones. The 
qualitative theory of the effect is presented. 
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Probing the mechanical response of single molecules 
to an external load has attracted considerable inter- 
est in recent years. In molecular failure experiments 
i|, i, S, i, H, 11 and experiments on protein unfolding 
7| a force grows linearly in time until a molecule breaks 
or changes its structure. The dynamic force spectroscopy 
0, Q delivering the spectrum of the rupture forces at dif- 
ferent loading rates gives deep insights into the internal 
dynamics of molecules 0, [ifl, [HI, In all these ex- 

periments polymers play an outstanding role, either as 
elastic couplers or as a subject of study 0, [H, 14 1. It 
was shown, that even the mechanical properties of pas- 
sive polymer spacers can affect the outcome of pulling 
experiments [5, 6, 14] . Recently a strong impact of the 

P'ymer size on its rupture dynamics has been reported 

The overall picture of the homopolymer mechanical 
failure is very similar to the one of the failure of macro- 
scopic fibers, as represented as a sequence of links. Under 
very slow ramping of the external load each link of the 
chain is subjected to the same force, and it is mostly 
the microscopic mechanism of a failure of a single link, 
which discriminates between the two situations. In a 
macroscopic case the single link failure is due to the pre- 
existing defects, so that the survival probability Wi{f) 
of a link under force / typically follows a power law. The 
probability that all the N links are intact is then given 

by 



WN{f) = W^[f{t)] 



N 



(1) 



and tends to a WeibuU law for fibers long enough. More- 
over, the longer fibers typically break at smaller forces, 
since the probability to find a weak link grows with the 
fibers' length, see e.g. [l6|. In a microscopic homoge- 
neous polymer chain, the single link breakdown is ther- 
mally activated; the breakdown probability at a given 
force at a given time follows an exponential law (as ob- 
tained from the Kramers theory of thermally activated 
barrier crossing), and therefore Wnif) tends to a Gum- 
bel distribution. However, the fact that the longer chains 
break more easily still holds: the larger is the number 



of links, the higher is the probability that one of them 
breaks [l3|. At higher loading rates the situation with 
the polymer chain gets more involved. As we proceed 
to show, a complex interplay between the temporal as- 
pects of thermally activated single link breakdown and 
the force redistribution along the chain leads to new fea- 
tures of the polymer failure problem. Thus, the behavior 
of the most probable force at breakdown as a function 
of the chain's length gets non-monotonic: the medium- 
length chains break more easily than the short and the 
long ones. The effect is not small and therefore is perti- 
nent to experimental observation. The breakdown force 
is reduced by 5 — 7% compared to the extremal cases of a 
single bond and of an infinitely long chain at the chosen 
values of the parameters, and the effect is expected to 
be the more pronounced the softer the bonds are since 
then the thermally activated breakdown happens earlier 
while the Rouse time, being a measure for the timescale 
of force propagation, becomes large. 

In what follows we first discuss our theoretical model 
and present the results of its numerical simulations. We 
then turn to the analytical description of the behavior 
observed, giving a full qualitative picture of the effect 
and then give its simple explanation. 

Like in [3] our model corresponds to a chain of N 
monomers interacting via the Morse potential 



(2) 



which parameterizes the interaction energy in terms of 
dissociation energy C / {2a) and stiffness Ca. This is a 
simple prototype of an intramolecular interaction poten- 
tial which offers fragmentation. Otherwise, the model 
is identical to the Rouse one [l3|: we disregard hydro- 
dynamical interactions and describe the interaction of 
the monomers with the heat bath via independent white 
noises. The constant loading enters through an addi- 
tional time dependent potential of the form 

L{q,t)^-qRt, (3) 

with loading rate R. The load is denoted by F(t) = 
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E"* — single bond rupture 
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Figure 1: Most probable rupture force jmax as a function of 
the chain length A''. The parameter values are c = 3.5nm/ ij.s, 
7? = 2 X 10"^nm^//is, r = 10"^nm/ps^, and a = lOnm"^ 
Error bars indicate the uncertainty due to the binning of nu- 
merical data. 

—dL{q, t)/dq and is applied at one end of the chain while 
the other end is fixed. 

At smaller loads the overall potential has two extrema, 

itxtri*) = ln(2)/a - In ± ^1 - F{t)/F)j /a, a mini- 
mum (jj^-tr corresponding to a metastable state of the 
pulled bond, and a maximum q^^i^ providing the activa- 
tion barrier. There exists a critical load Fc = F{tc) = 
C /4 for which the extrema merge at Qc = ln(2) /a and dis- 
appear. In the purely deterministic dynamics the Morse 
bond breaks exactly at = Fd R. Since the system is in 
contact to a heat bath at temperature T, its overdamped 
dynamics is described by a set of N coupled Langevin 
equations 

7(7i = -K{qi-qi-i) + K{qi+i-qi) + y/2k^Tj£,i + RtSi,N , 

(4) 

with K{q) = —dU{q)/dq, Gaussian white noise ^(t), the 
Boltzman constant k b , the friction coefficient 7 and go = 
0. We introduce c — C/7 ([c] = nm/ns), r = R/'f {[r] ~ 
nraj ^s^) and / = F (in the following / is referred to as 
force). The diffusion coefficient is denoted by £> = ksT 
{[D] = nm? / pLs). 

The set of coupled equations ^ is integrated by use 
of a Heun integration scheme. A chain is considered as 
broken as soon as one of the reaction coordinates qi+i — qi 
overcomes the activation barrier. Statistics stem from an 
ensemble of at least 10'^ simulation runs. 

In Fig. [1] we present the numerically obtained most 
probable rupture force fmax as a function of the chain 
length N for a fixed value of the loading rate r. For 
small chain lengths the numerically obtained most prob- 
able rupture force (symbols) follows the scaling relation 
given in Eq. ([7]) reaching a minimal value for an interme- 
diate value of A^. A further increase in the chain length 



Figure 2; Probability density of the position of the breakdown 
j in a chain of A*' = 100 bonds for different values of the 
loading rate r (given in the legend). The remaining parameter 
values are the same as in Fig. [T] 



results in an increase of the most probable rupture force 
and eventually in a saturation of the latter. A theoreti- 
cal description following the ansatz in [l^ (dashed line) 
gives in a good agreement the breakdown force for shorter 
chains as well as in the saturation regime, i.e. for very 
long ones. The curve obtained from Eq. ^ of the present 
contribution (solid line) gives the correct qualitative be- 
havior in the whole domain of the chains' lengths. 

The distribution of the position of the breakdown in a 
chain of A^ = 100 links is shown in Fig. [51 For a small 
loading rate, i.e., r = 10^^nm//is^ (solid line), there 
is only a very slight decrease of the rupture probability 
density along the chain. Virtually all bonds contribute 
equally to the rupture process. The situation changes 
drastically when passing to higher loading rates (dashed 
and dashed-dotted lines). The rupture probability den- 
sity decreases fast along the chain. For r = lO^^nm/ ^s^ 
only half of the chain contributes to the rupture process. 

In Fig. [3] we present the distribution of the position of 
breakdown in the chain for a fixed value of the loading 
rate and two different values of the length A^. One readily 
infers that the probability of a bond breakdown at a given 
site decays for a longer chain (solid line) faster with the 
distance from the pulled end than it does for a shorter one 
(dashed line). Thus, although the longer chains offers a 
larger number of possible breakdown sites, the number of 
bonds which contribute to the rupture process becomes 
smaller reaching a constant — loading rate dependent — 
value in the limit of a semiinfinitc system. 

In order to derive an analytical description of the 
chain rupture process we first recall the model describ- 
ing the single bond rupture. The probability Wi{t) that 
a single breakable bond remains intact can be expressed 
through the following kinetic equation 

0, Ei m 
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Figure 3: Probability density of the position of the breakdown 
j in a chain of A'^ bonds for one fixed value of the loading rate, 
r = W~^nm/ ns^ , and two different chain lengths (given in 
the legend). The remaining parameter values are the same as 
in Fig. [U 



dWiil 
rate 



dt 



-k(t)Wi{t), with k{t) being the Kramers 
Taking f(t) = rt we can rewrite the ki- 
netic equation in the form dWi{f)/df — '-k{f)/rWi{f). 
The measured probabiHty density function (PDF) for the 
rupture forces Pi(/) then is Pi{f) = -dWi{f)/df. 

Under the assumption that / is close to fc when bond 
rupture occurs, it is usual to expand the potential around 
the inflection point Qc up to the third order in deviations 
from [^, [13, 11 1 . We note that the breakdown prop- 
erties of the chain only depend on the behavior of the 
potential close to the point of critical load, which are 
universal U, 22, 2^. The Morse potential gives a conve- 
nient parametrization of the situation in terms of disso- 
ciation energy and stiffness. In the limit of small loading 
rates the most pro bable rupture force fmax follows the 

3 mi 



scaling relation 



fn 



In {v/r) 



(5) 



with V — ca^D/(87r) and w ~ c/{3aD). 

Passing to a chain of N breakable bonds we consider 
the rupture dynamics of different bonds to be indepen- 
dent. The probability that a bond i G 1 ... iV is intact is 
denoted by Wi{fi{t)). The pulling force at the chain end 
is / = /Ar(t). The equation for the survival probability 
of the chain reads: 



Wn (t) — exp ■ 



\n{Wi{f{x,t)))dx 



(6) 



If the typical rupture time tmax — fmax/r is much 
larger than the Rouse time r of the system, r = 
N'^ /{cair'^), (which is the case either for short chains 



Figure 4: Force profile in a Rouse chain with a coupling con- 
stant ca, t = The remaining parameter values are the 
same as in Fig. [1] 



or for small loading rates) each bond in the chain ex- 
periences the same force f{x,t) = f. The probability 
that a bond breaks in an interval [/, / + df] is = 
NWi{f)'^^^Pi{f) and it is given by the same expres- 
sion as Pi with V changed for Nv. This gives the scaling 
relation for the most probable rupture force [l5l | 



fmax fc 



In 



(7) 



The opposite situation is more involved. First, we have 
to calculate f{x,t). To do this, we note that the barrier 
crossing events are very rare: Most of the time the mo- 
tion of the monomers takes place close to the quadratic 
potential minima. Therefore, like in Ref. [l5| . the dy- 
namics of the chain can be approximated by a Rouse 
one. At difference to [l^ the boundary condition at the 
grafted end is explicitly taken into account. This leads to 
the dependence of the force profile on the chain's length 
as shown in Fig. [4l Linearizing the force profile close to 
the pulled end and inserting it into Eq. ([6|) we eventually, 
derive the PDF 



PiN,t)=.^WiN,t){^Ml), 

L3™3 5Ar(/)^ 

r(|a(/))-r0,a(/)(i + 5(iv,/))^ 

_ ^-dN{f)N ^_,(,)(i+5(Ar./))i 
rgNif) rgNif) 



with the probability function 



(8) 
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We introduced 
9N{f) = 
Scar 
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S{NJ) = gjv(/)iV/(/c - /), rfjv(/) = dgNif)/df and 
a(/) = w {1 — f / fc)'^^^ . r is the upper incomplete 
Gamma function. 

The analytical description given by Eq. ^ (solid line 
in Fig. [T])) agrees well with the outcome of the numerical 
simulations. The nonmonotonous behavior of the rupture 
forces is nicely reproduced and in the limit of short and 
long chains the theory agrees also qualitatively. Further- 
more the chain lengths which minimizes the most proba- 
ble rupture force coincide with the ones derived from the 
numerical simulations. Deviations can result from the 
harmonic approximation in the derivation of the force 
profile: the decay of the force profile along the chain for 
a soft Morse potential is expected to be more pronounced 
than in the Rouse chain, so that less bonds contribute to 
the rupture process. This might explain the shift of the 
theoretical curve to lower rupture forces compared to nu- 
merical data points. 

The overall behavior can therefore be explained as fol- 
lows. A short chain has only a few bonds that can break. 
Each of them feels practically the same force /, since 
the forces fi acting on the bonds decrease only slightly 
with their distance from the pulled end, see Fig. 21 The 
longer is the chain, the more breakable bonds are present, 
however each of them is subject to the tension which is 
smaller than / and decays with the chain's length. The 
interplay between the thermally activated rupture of a 
single bond and the force distribution along the chain 
generates a non-monotonous behavior of the typical rup- 
ture forces. 

Let us summarize our findings. Compared to the sin- 
gle bond breaking, the existence of the chain introduces 
new aspects into rupture dynamics, the most important 
being the delayed stress propagation along the chain. We 
show that the most probable rupture force decreases with 
the length of the chain as fmax oc — (ln(consi iV))^/'^ for 
short chains and saturates at the value depending on the 
loading rate for very long ones. In between it can exhibit 
a non monotonous behavior: the most probable rupture 



force attains its minimum for a certain intermediate chain 
length. These results are obtained via direct numerical 
simulations of a Rouse-like model (however with anhar- 
monic Morse interaction potential between the monomers 
instead of a harmonic one for a genuine Rouse chain) and 
are well reproduced by a generalization of an analytical 
approach put forward in our previous study The 
qualitative explanation of the effect involves a complex 
interplay between the force propagation into the chain 
and the extreme- value statistic underlying rupture. Since 
the effect is rather dependent on the different timescales 
in the system under study than on the specific parameter 
values it is pertinent to experimental observation. 
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